function y = func_delta2_approximated(p, m, CONSTS, plot_data)

    [q1, q2] = func_qk_from_p(p, CONSTS);
    k0    = CONSTS.k0;
    a     = CONSTS.a;
    Q1 = k0*a*q1;
    Q2 = k0*a*q2;

    y = besselj(m+1, Q2);

    if(plot_data)
        figure; plot(p, real(y), 'b-', p, imag(y), 'r-');
    end
    
end
